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Abstract 



Using factorization we determine the allowed domains of the ratios of form 
factors, x = A 2 (0)/Ai(0) and y = V(0)/Ai(0), from the experimentally 
measured ratio Rh = Y(Df — > <fip + )/Y(D+ — > <pTc + ) assuming three differ- 
ent scenarios for the g 2 -dependence of the form factors. We find that the 
allowed domains overlap with those obtained by using the experimentally 
measured ratio R S £ = Y(D+ — > <p£ + V()/T(D+ — > <pn + ) provided that the 
phenomenological parameter a^ is 1.23. Such a comparison presents a gen- 
uine test of factorization. We calculate the longitudinal polarization fraction, 
Yl/Y = Y(D+ — > (pLp\)/Y{D^ —* (f)p + ), in the three scenarios for the q 2 - 
dependence of the form factors and emphasize the importance of measuring 
Yl/Y. Finally we discuss the ^-distribution of the semileptonic decay and find 
that it is rather insensitive to the scenarios for the g 2 -dependence of the form 
factors, and unless very accurate data can be obtained it is unlikely to dis- 
criminate between the different scenarios. Useful information on the value of x 
might be obtained by the magnitude of the g 2 -distribution near q 2 = 0. How- 
ever the most precise information on x and y would come from the knowledge 
of the longitudinal and left-right transverse polarizations of the final vector 
mesons in hadronic and/or semileptonic decays. 
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Since the Cabibbo-Kobayashi-Maskawa(CKM)-favoured decays of the strange-charmed 
meson Df involve, at the quark level, ss pair, the vector meson <ft emerges naturally as 
one of the decay products. Among its decay modes are the three most important ones, 
Df — > 4>n + , D+ — > cj)p + and Df — > <f)£ + i>i(£ = e, p) measured by several groups Q, and 
the topic of this paper. 

The analysis of the experimental data for these modes involves three axial form factors, 
Aq, A\ and A%, and a vector form factor V, in the notation of Bauer, Stech and Wirbel 
(BSW) [0]. In spite of this fact which renders the analysis of experimental data involving 
the meson more complicated than the corresponding one with pseudoscalar mesons 
such as 7] or r/, it turns out that the polarization of the <fi meson may provide crucial 
information that can be used to test the factorization assumption which is one of the 
basic ingredients in the computation of hadronic decay rates. 

The fact that the n meson is light allows us to ignore the variation of the hadronic form 
factor A between q 2 = m 2 , where it enters the description of Df — > 07r + , and q 2 = 0. 
This then allows us to eliminate A s< ^(0) in favour of A x (0) and A 2 "^(0) through 0, 

2m^(0) = (m Ds + m^Af'^O) - (m Ds - m^)A^(O) (1) 

Since the other decays, Df — > (pp + and Df — > (p£ + U£, involve the form factors A Y , 
A 2 and V Da ^, it proves convenient to factor out A x (0) which then cancels in the ratio 
of the rates, and define the following ratios of the form factors at q 2 = 0, 

,4^(0) _ V D ^(0) 



X ^ >D s ct>,n\ y ~ aD s S,^ (2) 



In our analysis we do not adopt any of the several models of hadronic form factors for 
D+ — > transition available in the literature to compute decay rates as such a procedure 
would not tell whether the factorization assumption or the model-form factors were being 
tested in the comparison between theory and experiment. Instead, we simply consider 
different scenarios for the g 2 -variation of the form factors ( and do not make any assump- 
tions about their absolute normalizations ) and show that the measured ratio of the two 
hadronic rates, 

R h = T(Df -> <P P + )/T(DJ -> 07T+) (3) 

constrains the parameters x and y to an allowed domain of the x, y plane ( limited by 
x and y non-negative a suggested by all phenomenological models ). Such a domain can 
then be compared with the values of x and y already measured in the semileptonic modes 
P], [51 01]. We remind the reader that semileptonic data are analyzed with monopole q 2 - 
dependence for all form factors which is just one of the scenarios we discuss. A procedure 
such as ours that eschews reliance on a particular set of model-form factors produces not 
only a more bias-free test of factorization but also one that is independent of the BSW 
phenomenological parameter a\ @. 

Further, since the <ft meson is longitudinally polarized in D^ —>■ 4>tt + , the ratio T(Df — > 
<f)7T + )/r(D+ — * <pp + ) behaves physically like the longitudinal polarization fraction Yl/Y 
in Df — > 4>p + , where 



s 



r Y(D; 



(4) 



Therefore, there must be a link between these two quantities Rh and Yl/Y which we 
explicitly demonstrate. Indeed, an upper and a lower limit for Y^/Y can be calculated for 
the allowed x, y domain. A direct measurement of this quantity would provide a further 
test of factorization. 

We now recall, for completeness, the various expressions for decay rates used in this 
paper though they can be found in published literature || ||, 0]. We start with the 
semileptonic decay Df — > 0£ + z^ in the zero lepton mass limit. Let us denote by q 2 the 
(effective mass) 2 of the ^-system and introduce the following dimensionless quantities, 

r = — — t = —5- (5) 

m Da m z Ds 

with < t 2 < (1 - rf. 

In the Df rest system, the meson momentum K(t 2 ) is given by, 

K(t 2 ) = ^f k(t 2 ) (6) 

where k(f) = [(1 + r 2 - t 2 ) 2 - 4r 2 ] 1/2 - 

The £ 2 -distribution for semileptonic decay is the sum of the contributions from the 
three polarization states, 

|j (Df - <pv t ) = ^|^ |V CS | 2 k{t 2 ) {H LL (t 2 ) + H ++ {t 2 ) + H^{t 2 )} (7) 



where 

i-'^\ _ l' ~ r ' V- It -1 ..2 .,2, i/i--;.,.:?, fe (^ ) AD a (f>/j2\\2 



"^na-^-^^v)-^ 



^LL(t 2 ) = [^] 2 |(l-r 2 -tVf*V)- 



tf ±± (t 2 ) = (l + r) 2 t 2 lA^^TjPp^V^it 2 )] 2 (8) 

For the hadronic modes where only the spectator diagram contributes in the factor- 
ization approximation, we obtain 

Y(Df -> 0vr+) = %^ IK.IW a 2 (^ ) 2 MO) ^ll(O) (9) 

OZ7T TYlDs 

where for t = 0, 

MO) H LL (0) = (1 - r 2 ) 3 [i±I ] 2 {1 - i^I a;} 2 |A?'*(0)| 2 (10) 

For the 0p + mode, first treating the p meson as a zero-width resonance, we get 

Y(Df -> #+) = %^ IVcHV^I 3 a? (^ ) 2 fc© {tf LL (t 2 ) + H ++ {t 2 ) + H&)} 



(11) 
where t p = m p /mD s - A comparison of (9) and (11) provides us with the precise connection 
between Y(D+ — > </>7r + ) and r(_D+ — ► 4>lPl)i ^ ne ra te for the p meson in the longitudinally 
polarized state. 



To take the finite p-width into account one has to smear the rate given in (11) over t 2 
with a Breit-Wigner measure 0, 

*(# ff «© — \ j^ dt \ t ._lp +ipl H^H^) (12) 

where 7 P = T p /rriD a , T p being the p- width. 

The analysis of the experimental data requires certain assumptions for the g 2 -dependence 
of the hadronic form factors. In our analysis we assume the vector form factor, V(q 2 ), to 
have a monopole g 2 -dependence with a pole mass mi- = 2.11 GeV M, As for the axial 
form factors, Ai(q 2 ) and A 2 (q 2 ), we investigate the following three scenarios : 

Scenario I : Both Ai(q 2 ) and A 2 (q 2 ) assumed to have monopole g 2 -dependence with 
pole mass ra\+ = 2.53 GeV as in BSW @. 

Scenario II : Both Ai(q 2 ) and A 2 (q 2 ) assumed to be independent of q 2 in the region of 
interest. 

Scenario III : A\(q 2 ) assumed to be decreasing with q 2 , 

A^(q 2 ) = Af'*(0)[l - (3 - q2 QeV)2 } (13) 

while A 2 (q 2 ) is assumed to have a monopole dependence with a pole mass of 3.5 GeV. 
This latter value was chosen just for illustration of a gently decreasing and increasing 
Ai(q 2 ) and A 2 (q 2 ) respectively. Qualitatively, Scenarios II and III were inspired by the 
QCD sum rule calculations of Ref. [§ for D — > K* transition form factors. 
Using the measured hadronic ratio, 

* = n%Z%?) = 1M±0 - x %£ (14) 

we calculated the allowed domains in the x,y plane for each of the above scenarios. For 
T(D+ — > 07r + ) we used expressions (9) and (10), and for T(D+ — > <fip + ) the t 2 -smeared 
form of (11) in accordance with the prescription of Eq.(12). 

In figs. 1, 2(a) and 2(b) we have plotted the boundary curves for scenarios I, II and 
III respectively. All positive values of x and y are allowed up to the boundary curve. 
Experiments 0, [| on semileptonic decays, Df — > 0£ + z/£, have determined x and y using 
the assumptions of scenario I. We have plotted their points in Fig. 1 too. We note that 
the point from E-653 || is excluded by the constraint set by Rh while the point from 
E-687 [Q has a small overlap with the domain determined from R^. In our analysis the 
two errors in (14) were added in quadrature and we used the leptonic decay constants 
f w+ = 131.7 MeV and f p + = 212 MeV in Eqs. (9) and (11). 

Technically the curves for fixed Rh are hyperbolae in the x, y plane centered on the x 
axis whose principal axes are the x axis and an axis parallel to the y one. Of course only 
the part of the curve located on the left branch of the hyperbola has a physical meaning; 
the second branch requires too large a value of x, x > 5.0. 

The general equation for these hyperbolae is, 



I iAj iA-J f> 



a L 



y_ 

b 2 



^ =1. 



(15) 



We have listed the numerical values of x c , a and b for the three models in Table 1 where 



we have also tabulated the upper limits of x and y, 
and, of course, x max — x c - a. 



x i 



for y = and y max for x = 0, 



SCENARIO 


x c 


a 


b 


Xmax 


Vmax 


I 


2.81 


2.37 


6.14 


0.44 


3.92 


II 


2.87 


2.13 


5.57 


0.74 


5.01 


III 


2.89 


2.00 


5.20 


0.89 


5.45 



Table 1. 
Parameters of the hyperbolae and x max and y ma x of the allowed x, y domain. 

The full curves in Figs. 1, 2(a) and 2(b) use the smeared rate for T(Df — »■ <fip + ) in 
Eq.(12) while the dotted curves use the p + zero-width approximation in Eq.(ll). We 
remark that for scenario I the allowed x, y domain for the p + zero-width approximation 
is empty. We emphasise that the corrections due to the non-zero p meson width play an 
important role in all scenarios. 

Having established the range of x and y allowed by Rh, we can calculate using (11) 
and (12), the maximum and the minimum values of the longitudinal polarization fraction. 
We find the following limits; 



scenario I : 
scenario II : 
scenario III : 



0.43 < -f? < 0.55 

■p 
0.36 < -jr < 0- 55 

■p 
0.33 < -f? < 0.55 



(16) 



The upper limit occurs for x = y = and is insensitive to the scenario adopted, while 
the lower limit depending on the maximum value of y, which occurs at x — 0, is rather 
scenario-dependent. In passing, we note that for a zero- width p-meson the maximum 
value of T L /T is 0.51. 

It is evident that an experimental measurement of the vector meson polarization in 
the decay D+ — > <pp + could provide important information on the parameters x and y. 
Further, a large experimental value of T L /Y exceeding unambiguously the upper bound 
given in (16) will be a clear indication of the failure of the factorization approximation. 
Equally, a measurement of Tl/T below 0.33 would rule out all the three scenarios assum- 
ing, of course, the validity of factorization approximation. In short, Tl/T could be used 
to discriminate between different scenarios. 

In Figs. 1, 2(a) and 2(b) we have also plotted constant Tl/T curves (dashed-dotted) 
for three values, 0.3, 0.4 and 0.5 in order to exhibit the variation of Tl/T in the x, y 
plane. As previously, these curves are also hyperbolae of the type given by Eq.(15). 



We now turn to a discussion of the semileptonic decays Df — > (f)d, + v^. Experimentally, 
the following ratio has been measured |1|, £|, [1(| [11], [L2| 

R, ^ r ^ Tn = °- 54±0 - 10 < 17 > 

r(-D+ -► 07T + ) 

Using (7) and (8) we can calculate T(Df — > 4>(i + U() for each of our three scenarios 
for the g 2 -dependence of the form factors. Then using (9) and (10) we can construct a 
theoretically ax-independent ratio, 



T(D+ - <j>Fv t ) _ 2 



T(D+ - 0ttH 



i -Rs 



If ai were a well-determined parameter we could use (18) with R S £ given by (17) to 
determine the allowed x, y domain from the ratio R S £. It turns out that the limits of 
the allowed domains are rather sensitive to the value of oi used. If we treat a± as a free 
parameter we find that for Oi ~ 1.23 the allowed domains obtained from Rh and R s e 



almost overlap. This is quite remarkable as the expected value of a\ is 1.26 [I3|. This we 



interpret as an evidence in favour of the factorization hypothesis in the hadronic decays. 

We now address the question : What can we learn from the g 2 -distribution, dT(D+ —+ 
(f>£ + b>e) J dq 2 , of the semileptonic decays ? 

For 0~ — > 0~£i>f decay the answer is simple - we determine the g 2 -dependence of 
the form factor Fi(q 2 ) as there is only one form factor involved. In the present case of 
CT — * l _ £z/ decay the answer is not so clear as there are four terms in the g 2 -distribution 
associated with A 2 (g 2 ), Ai(q 2 ) ■ A 2 (q 2 ) , A 2 ,^ 2 ) and V 2 (q 2 ) and the only way to disentan- 
gle these four contributions, and then to measure the g 2 -dependence of the various form 
factors, is to determine simultaneously the longitudinal and left-right transverse polariza- 
tions of the final vector meson. We are far from achieving this in the case of D+ -^ <fi(i + vt. 
decay in the near future. 

Unfortunately, experimental measurement of the true g 2 -distribution of the semilep- 
tonic decays has not yet been possible because the momentum of the D^ decaying in flight 
is not known |14| which, in turn, implies that the q 2 for an event is not accurately known. 
However, two experiments, E-653 and E-687 have determined the parameters x and 
y, assuming scenario I for the g 2 -dependence of the form factors, from the semileptonic 
decay Df — > ipl + vi. In the following we exploit this knowledge in absence of information 
on the g 2 -distribution dT /dq 2 . 

Let us introduce a unit normalized dimensionless quantity X(t 2 ), 

x(t 2 ) = i ^ rp+ - <ptu,) (19) 

which depends on x and y, and on the type of g 2 -variation of the hadronic form factors. 
We first construct an "experimental" X(t 2 ), denoted by X(t 2 ) in the following, using the 
experimental values of x and y obtained from a triple angular fit of the data 0, |j] in the 
scenario I. We then compute the same quantity, X(t 2 ), for the values of x and y inside the 
allowed domains determined by us from R^ and shown in Figs. 1, 2(a) and 2(b). Finally 
we compare the X(t 2 ) so constructed with X(t 2 ). 



From Fig. 1 we note that the E-653 point lies well outside the allowed x, y domain 
while the E-687 point has a small overlap with this domain. Consequently we use only 
the latter data point in constructing X(t 2 ) and we adopt [|J, 

x = 1.1 ±0.81, y = 1.8 ±0.92 (20) 

using scenario I, as was also the case in j||. 

The resulting X(t 2 ) is shown in Fig. 3(a) in the form of nine curves corresponding to 
the central and one standard deviation extreme values of x and y. These nine curves break 
up naturally into three groups of three curves each, with each group being characterized 
by the value of x ; within each group the splitting due to different values of y is a much 
finer effect. Qualitatively at t 2 = 0, the group with the smallest x has the largest intercept 
with the vertical axis and within each group the curve with the smallest y has the largest 
intercept. 

Our calculated X(t 2 ) for scenarios I, II and III are shown in Figs. 3(b), 4(a) and 4(b) 
respectively in the form of three curves corresponding to the three extreme values of x and 
y for the allowed domains of Figs. 1, 2(a) and 2(b), (x, y) = (0, 0), (0, y ma x), (x max , 0). All 
other choices of x and y in the allowed domain yield X(t 2 ) which lies within the envelope 
of the three curves shown in Figs. 3(b), 4(a) and 4(b). The t 2 -distribution in these figures 
are then our prediction for X(t 2 ) for x and y obtained from the hadronic ratio R%. We 
observe that they depend strongly on x and very little on the q 2 scenarios. Therefore, we 
conclude that X(t 2 ) cannot distinguish between different q 2 scenarios unless very precise 
measurements of the ^-distribution can be made which probably will not be the case in 
the near future. However, the determination of X(t 2 ) near t 2 = appears to be very 
useful for the determination of the parameter x. 

We now compare our curves X(t 2 ) with the "experimental" ones, X(t 2 ), of Fig. 3(a). 
For scenario I, we find some overlap between the curves of Fig. 3(a) and 3(b) provided 
that x is small. In particular the group of curves in Fig. 3(a) with x = 0.29 overlap well 
with the curve in Fig. 3(b) which uses x = x max and y = 0. The reason is obvious from 
Fig. 1 which shows that the points, x = 0.29 and y between 1 and 2 lie inside the allowed 
domain of x and y obtained from R^. We remind the reader that X(t 2 ) depends only 
weakly on y. 

Such an overlap appears a little more problematic in scenarios II and III especially for 
large values of t 2 but cannot be excluded. For this reason we suggest to experimentalists 
to fit their triple angular distributions with different scenarios for the g 2 -dependence of 
the hadronic form factors for extraction of x and y. It would be interesting to see if the 
corresponding extracted values of x and y for scenarios II and III would be consistent 
with the allowed domains of x and y for these scenarios extracted from R^ and shown 
respectively in Figs. 2(a) and 2(b). 

In conclusion, assuming factorization and using R h (see(3) and (14) for definition), 
appropriately smeared over the p-width, we found the allowed domains in the x, y plane 
using three different scenarios for the g 2 -dependence of the form factors. These allowed 
domains plotted in Figs 1, 2(a) and 2(b) show sensitivity to the choice of the scenario. 
There is an overlap, though small, with the E-687 || values of x and y from semileptonic 
decay T(Df —>■ <ft£ + i>£) and the x, y allowed domain calculated by us from R^ in the 
relevant scenario - scenario I ( See Fig. 1). This would suggest that the hadronic data 



from H and the semileptonic data from || are consistent with the factorization hypothesis. 
E-653 || determination of x and y from semileptonic decay, on the other hand, does not 
support factorization assumption. 

We have also determined the allowed domains of x and y from R s i (see (17) for 
definition) in the three scenarios. Such a determination is ax-dependent. We found that 
the two domains determined from Rh and R s e overlap almost completely for oi « 1.23 in 
remarkable agreement with the theoretically computed value of a\ = 1.26 [[H|]. This fact 
supports the validity of the factorization assumption. 

We have plotted the contours of constant Y L /Y ( see (4) for definition) in the x, y plane 
and calculated the maximum and the minimum values of Yl/Y for the three scenarios of 
g 2 -dependence of the form factors considered by us. It is remarkable that the maximum 
value of Yl/Y = 0.55 is independent of the scenario for the g 2 -dependence of the form 
fators, while the minimum value is scenario-dependent. Of course, a measurement of the 
longitudinal and left-right transverse polarizations will determine both x and y in a chosen 
scenario for the g 2 -dependence of the form factors. 

A measurement of Yl/Y could be important in testing factorization. For example, 
should the measured value of Yl/Y exceed 0.55 unambiguously, it would signal a break 
down of the factorization hypothesis. A value below 0.33 would rule out all of our three 
scenarios for the g 2 -dependence of the form factors if factorization were to hold. 

Finally, we suggest that experimentalists analyze their semileptonic data using for q 2 - 
dependence of the form factor forms other than monopoles to extract the value of x and 

y- 
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Figure captions 



1. Fig. 1 : The shaded region bounded by the solid line is the allowed domain of x, 
y from R^ within one standard deviation. The dashed-dot lines correspond to fixed 
values of Tl/T as indicated. All these curves are calculated with the scenario I for 
the q 2 dependence of form factors. 



2. Fig. 2(a) and 2(b) : Same as Fig. 1, for scenarios II and III respectively. The 
dotted lines replace the solid lines in the zero width approximation of the p meson. 



3. Fig. 3 : 

(a) "Experimental" X(t 2 ) as explained in the text. 

(b) Calculated X(t 2 ) in the scenario I. 

4. Fig. 4 : 

(a) Calculated X(t 2 ) in the scenario II. 

(b) Calculated X(t 2 ) in the scenario III. 
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